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Abstract 

We study the effect of coherent propagation of two interacting particles in 

an effective 2-3-d disordered potential. Our numerical data demonstrate that 

in dimension d > 2, interaction can lead to two-particles derealization below 

one-particle derealization border. We also find that the distance between 

the two delocalized particles (pair size) grows logarithmically with time. As 

a result pair propagation is subdiffusive. 

PACS numbers: 71.55.Jv, 72.10.Bg, 05.45.+b 
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I. INTRODUCTION 



The question of interacting particles in a random potential has recently got a great deal 
of attention (see for example 0). Indeed this problem is important for the understanding 
of conduction of electrons in metals and disordered systems. It is also very interesting 
from the theoretical view point since it allows to understand the effects of interaction on 
Anderson localization. It is a common belief that in one-dimensional (Id) systems near 
the ground state, a repulsive interaction between particles leads to a stronger localization 
if compared with non interacting case 0. Even if more complicate, the 2-dimensional (2d) 
problem is assumed to be localized, while in the 3-dimensional (3d) case derealization can 
take place in the presence of interaction 0]. However this problem is rather difficult for 
analytical, experimental and numerical investigations and therefore it is quite far from its 
final resolution. 

The complicated nature of the above problem can be illustrated by the example of only 
two interacting particles (TIP) in a random potential. Indeed in this case, contrary to the 
common lore, even repulsive particles can create an effective pair which is able to propagate 
on a distance l c much larger than its own size of the order of one particle localization length 
l\ 0. From one side interference effects for non-interacting particles force the particles to 
stay together at a distance ~ l\ even in the repulsive case. But the relative motion between 
the two particles leads to the destruction of such interference and allows their coherent 
propagation on a distance l c ^> l\. More explicitly, according to |§, in the quasi Id case 
with M transverse channels one has 

l tr hM ^ (L1) 

where U is the strength of on site interaction and V is the one particle hopping matrix 
element. Here the inter-site distance is a = 1 and the wave vector kp ~ 1. Since l\ oc M, 
then l c oc M 3 . In 2d the effective number of channels due to 2d localization is M ~ l\ so 
that l c ~ l\ and localization is preserved (here and everywhere li represents the one-particle 
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localization length in any dimension). The sharp increase of l c with the number of transverse 
channels M leads to a straightforward possibility of derealization for a pair of particles in 3d 
while one particle remains localized |§-|6|]. In this sense the 3d case is much more interesting 
due to the possibility of derealization and we qualitatively discuss it below. 

One of the interesting features of pair derealization in 3d is that it is not due to a 
simple shift of the mobility edge produced by the interaction (such possibility is indeed 
not so interesting). In fact, it is possible to consider a system in which all one-particle 
eigenstates are localized for all energies. This can be for example the 3d Lloyd model 
with diagonal disorder _E„ 1>n2in3 = tan0 ni)n2jn3 and hopping V on a cubic lattice, where 
4>m,m,nz are random phases homogeneously distributed in the interval [0,7r]. In this model 
all one-particle eigenstates are localized for V < V c ~ 0.2. However, two repulsive particles 
with on-site interaction U can create a coupled state which is delocalized and propagates 
through the lattice. Indeed, a pair "feels" only smoothed potential || that corresponds to 
an effective renormalization of the hopping matrix element V e ff which is strongly enhanced 
due to interaction and becomes larger than V c . Of course, the enhancement takes place only 
for sufficiently large one-particle localization length / x 3> 1. Therefore, the hopping V < V c 
should be not far from V c although there is, a priori, no requirement for V to be very close 
(parametrically) to V c . 

The two particles derealization due to interaction takes place only for states in which 
particles are on a distance R < l\ from each other while for R 3> l\ eigenstates are localized. 
Such kind of situation is quite unusual since it means that the absolutely continuous spec- 
trum of the Schrodinger operator, corresponding to the delocalized pair, is embedded into 
the pure point spectrum of localized almost noninteracting particles states. For R^$> l\ the 
interaction between the two particles is exponentially small and this implies a very small 
coupling between the states corresponding to these two kinds of spectra. However, due to 
the quasi degeneracy of levels, even a small coupling can lead to important modifications 
of the above picture as it was discussed in |J. Therefore, a direct numerical investigation 
of interaction-assisted derealization in 3d is highly desirable. While the recent theoretical 
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arguments and numerical simulations in quasi-ld case - definitly demostrate the ex- 
istence of enhancement for l c no numerical simulations have been done in 3d case. Indeed, 
in 3d basis grows as iV 6 , where iV is the number of Id unperturbed one-particle states, and 
that leads to heavy numerical problems. 

A similar type of interaction-assisted derealization can be also realized in the kicked 
rotator model (KRM) JlO| in 3d. In this case the unitary evolution operator takes the place 
of Schrodinger operator and eigenenergies are replaced by quasi-energies. The advantage 
of such models is due to the independence of localization length on quasi-energy so that 
all one-particle states in 3d are localized for V < V c and delocalized for V > V c . Even if 
very efficient, numerical simulations for KRM in 3d become very difficult; for two particles 
situation becomes even worse due to N 6 basis growth. 

One of the ways to overcome these numerical difficulties is the following. For Id KRM 
the number of dimensions can be effectively modelled by introducing a frequency modulation 
of the perturbation parameter [Tl||l^]. The case with v incommensurate frequencies in the 
kick modulation corresponds to an effective solid state model with dimension d — v + For 
v = 2, the effective dimension d = 3 and Anderson transition can be efficiently investigated 
\T2\. A similar approach can be done for two interacting particles and it allows to gain a 



factor N A in numerical simulations. 

In this paper we investigate the model of two interacting kicked rotators (KR) studied 
in 0, with frequency modulation and v = 2,3. The quantum dynamics is described by 
the evolution operator : 

S 2 = exp{-i[H (n) + H (n') + U8 n>n ,\} ^ 
xexp{-i[V(9,t)+V(9',t)}} 

with = —id/d9 ( -'\ Here H (n) is a random function of n in the interval [0, 2n] 
and it describes the unperturbed spectrum of rotational phases. The perturbation V 
gives the coupling between the unperturbed levels and has the form V{9, i) = k(l + 
e cos 6*i cos 6 2 cos 9 3 ) cos 9 with #1,2,3 = ^1,2,3 t. In the case of two modulational frequen- 
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cies [y = 2,uj 3 = 0), as in ||12|| , we choose frequencies ui^ to be incommensurate with each 
other and with the frequency 2n of the kicks. Following [[L^] we take uji = 27rA _1 , u>2 = 2tt\~ 2 
with A = 1.3247... the real root of the cubic equation x 3 — x — 1 = 0. For v = 3 we used 
the same uji^ and cu 3 = 2tx j \[2. We also studied another case of functional dependence 
of V(9) analogous to ]I2| and corresponding to the Lloyd model. All computations have 
been done for symmetric configurations. According to the theoretical arguments || and 
numerical simulations || the antisymmetric configurations corresponding to fermions with 
nearby site interaction should show a similar type of behaviour. 

The paper is constructed as follows. In section II we discuss the model and present the 
main results for v = 2. The case of v = 3 is discussed in section III. The kicked rotator model 
corresponding to the 3d Lloyd model is studied in section IV. Conclusions and discussions 
of results are presented in section V. 

II. THE KRM MODEL WITH TWO FREQUENCIES 

Before to discuss the effects of interaction let us first discuss the noninteracting case 
U = 0. Here the evolution operator can be presented as a product of two operators Si 
describing the independent propagation of each particle: 

Si = exp(-iH Q (h)) exp(-iV(9, t)) (2.1) 

Since V depends on time in a quasiperiodic way with V(9, t) = V(9, 9 X , 6*2) and #12 = uji^t, 
one can go to the extended phase space |lTJ, [12| with effective dimension d = 3. In this 



space the operator is independent on time and has the form 

Si =exp(-iH 1 (h 1 h 1 ,n 2 ))ex P {-iV(e,e li e 2 )) (2.2) 

with Hi(n, m, n?) = Hq(ti) + u)in\ + ^2^2- Due to linearity in ni^ the transformation from 
( p.ip to (|2.2|) is exact. However, the numerical simulations of ( |2.1| ) are iV 2 times more 
effective than for (]2.2). 
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The system ( |2.2|) corresponds to an effective 3d model. Numerical simulations in [12 



showed that the variation of coupling amplitude V gives the transition from localized to 



diffusive regime as in usual Anderson transition in 3d. In |L2| the form of the kick V had 
been chosen as 

V(d,9i,9 2 ) = -2tan- 1 [2A;(cos^ + cos^i + cos^2) - E\ (2.3) 

In this case after a mapping similar to the one used in [|H| the equation for eigenfunction 
with a quasi-energy /i can be presented in a usual solid-state form: 

T n u n + kJ2 r u n - r = Eu n (2.4) 

where the sum is taken only over nearby sites and T n = tan((.£Zi(n, ni, n 2 ) — aO/2), 
n = (n, ni,n 2 ). For random phases under tangent the diagonal disorder is distributed 
in Lorentzian way and the model becomes equivalent to the 3d Lloyd model. While we also 
investigated the kick form ( |2.3j ) (see section IV) our main results have been obtained for 

V(0, 0i, 9 2 ) = k cos 6(1 + e cos 6 l cos 2 ) ( 2 - 5 ) 



in the case of two frequencies v = 2. According to [TJ| in this case the equation for eigen 



functions can be also reduced to an effective solid-state Hamiltonian which, however, has 



a bit more complicated form than fl2.4j ). We choose (|2.5| ) since it was numerically more 
efficient than (|2.3|) . To decrease the number of parameters we always kept e = 0.75. 

The one-particle transition as a function of coupling (hopping) parameter k in (|2.5| ) is 
presented in Fig.l. Similar to |12j the localization length li is determined from the stationary 
probability distribution over unperturbed levels \ip n \ 2 ~ exp(— 2|n|//i) while the diffusion 
rate is extracted from the gaussian form of the probability distribution lnW^ ~ —n 2 /2Dt 
with D =< n 2 > jt. According to Fig.l the transition takes place at the critical hopping 
value k cr w 1.8. Below k cr all quasi-energy states are localized. The independence of 
transition point from quasi-energy is one of useful properties of KR models. 

Our main aim was the investigation of TIP effects well below the transition point k cr . 



As in H we characterized the dymamics ( |1.2|) by the second moments along the diagonal 

6 



line n = n' : u+(t) = ((\n\ + |n'|) 2 ) i /4 and across it cr_(t) = — \n'\) 2 ) t . We also 
computed the total probability distribution along and across this diagonal |J: P±(n±) with 
n± = \n ±n'\/2 1 / 2 . The typical case is presented in Figs. 2,3. These pictures definitely show 
the appearence of pair propagation even if the interaction is neither attractive nor repulsive. 
Indeed, the pair size is much less than the distance on which two particles are propagating 
together. If we fit the probability distibution in Fig. 3 as P± ~ exp(— 2n±/l ± ) then we can 
see that the ratio ~ 25 (l c ~ / + ~ 95) is quite large. It is interesting to note that 

P + {n + ) at different moments of time is closer to an exponential (lnP + ~ n + ) than to a 
gaussian (lnP + ~ n + 2 )- The spreading along the lattice leads only to growth of l + with 
time but the shape of distribution does not corresponds to a diffusive process. 

Another interesting feature of Fig. 2 is the slow decrease of the rate of a + growth and 
the slow growth of o~_. To check if the growth of a + is completely suppressed with time 
we analysed its dependence on t for different values of interaction U (Fig.4, probability 
distribution is shown in Fig.5). For U < 0.5 the growth of a + is completely suppressed 
while for U = 1 the complete suppression is a bit less evident. To understand in a better 
way the case U = 1 we can look on the dependence of the number of effectively excited 
levels AiV on time. To estimate AiV we should rewrite ( |1.2j ) in the extended basis where 
the evolution operator has the form: 

S 2 = exp(-t[H (n) + H (n') + uj.n, + oo 2 n 2 + U5 n y}) exp(-i[V(9, 6 U 6 2 ) + V(9', 6 Xl 6 2 )\) 

(2.6) 

Since Ani^ ~ An^ « An + the number of excited levels can be estimated as AiV « 
An + An_AniAn2 ~ Following the standard estimate based on the uncertainty 



relation |TJ| a derealization can take place only if AiV grows faster than the first power of 
t. As our numerical data show ( see Fig. 6) the ratio W = AN/t remains approximately 
constant or is even slightly decreasing in time. This indicates that our case is similar to 
localization in 2c? where this ratio also remains constant for very long time. The reason why 
below k cr the situation is similar to 2c? can be understand in the following way. According 



to ( |2.6|) the total dimension is 4 and we have there 2 particles. Therefore, we can argue 
that the dimension per particle is 2 and that below the 3d derealization border k cr our 
system effectively represents two particles in an effective dimension d e ff = 2. However, two 
particles in 2d are always localized but the localization length can be exponentially large. 
The dependence of cr + at fixed moment of time on k is presented in Fig. 7 and indeed, it 
demonstrates a sharp increase of cr+ with l\ approaching k cr . Therefore, we conclude that for 
k < k cr our model effectively represents TIP in 2d. Below k cr the pair created by interaction 
remains localized but the localization length l c grows exponentially with l\. To see effects 
of interaction in d e ff > 2 we should study the system with three modulational frequencies 
v = 3. But before to analyse the case v = 3 we would like to discuss the behaviour of cr_. 

Indeed, Figs. 2,4 clearly demonstrate a slow growth of cr_ with time which means the 
increase of the size of the pair k. The results presented in Fig. 8 show that k m a^ 1 ! 2 grows 
logarithmically with time as k ~ C^lnt where Cl is some time independent factor being 
Cl ~ 0.8(£7 = 2) and Cl ~ O.Q(U = 1). Of course, this logarithmic growth should terminate 
after the complete localization in a + but this time scale t c is very large and for t < t c we have 
clear logarithmic growth of k. As discussed in || we attribute this growth to the fact that 
propagating in a random potential the pair is affected by some effective noise which leads 
to a slow separation of two particles. Indeed, the matrix elements of interaction U- decay 
exponentially fast with the growth of the pair size n_ = k according to a rough estimate 
[/_ ~ U s exp(— \n-\/li) with U s ~ V ' jl\^ 2 '. These small but finite matrix elements lead to 
the growth of the pair size k with slow diffusion rate D_ oc U s 2 exp(— 2|n_|// 1 ). According 
to the relation n 2 /t ~ D_ the pair size grows as k ~ l\ lnt/2 which is in agreement with 
data of Fig. 7. More detailed numerical simulations are required to verify the dependence 
Cl ~ h- 
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III. THE KRM MODEL WITH THREE FREQUENCIES 



According to the above discussion the suppression of diffusive growth of cr + can be 
explained by two factors. The first one is that the effective dimension is d e ff = 2 and 
localization always takes place in 2d. Another reason is the slow logarithmic growth of pair 
size. To separate these two effects we investigated the dynamics of TIP in the KRM with 
three modulational frequencies v = 3. In the extended phase space the evolution operator 
has the form 

S 2 = exp(-i[H (n) + H (n') + uj 2 n 2 + u 3 h 3 + U5 n>n/ ]) 

exp(-i[v(e, e u 2 , e 3 ) + v(e', e ± , e 2 , e 3 )}) 

with 

V(6,6 U 6 2 ) = k cos 6(1 + e cos 9 X cos 6 2 cos 9 3 ) (3.2) 

For U = we have one particle in 4d and transition to derealization takes place above a 
critical value of perturbation parameter k cr . According to our numerical data k cr m 1.15 for 
e = 0.9 (Fig. 9). Below k cr all eigenstates are exponentially localized. For (7^0 the total 
dimension in ( |3.1| ) is 5 and since we have 2 particles the effective dimension per particle is 
d e ff = 5/2. Since d e ff > 2 the first argument given above becomes not relevant and we 
expect TIP derealization below k cr . Let us to note that above k cr the above TIP problem 
becomes not interesting since even without interaction the particles spread along the lattice 
and the interaction between them does not affect significantly their dynamics. 

The numerical simulations of TIP for (|3.1| ) - (|3.2| ) in one-particle localized phase k < k cr 
demonstrate strong enhancement of two particles propagation. A typical case is presented 
in Figs. 10,11. According to these data the growth of a + is unlimited and TIP dereal- 
ization takes place below k cr . The analysis of cr_(t) shows that pair size k ~ cr^ 1 / 2 grows 
logarithmically with time similar to the case with v — 2. We think that this slow growth 
of k is responsable for a slow decrease of the pair diffusion rate D + = a + /t with time. Due 
to the increase of k the probability to have a distance between particles of the order of l\ 
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decreases as l/n{t) ~ 2/(/ilnt) and therefore, we expect that the diffusion rate of the pair 
will decrease with time as D + ~ D e f/\r^t. Here \x = 1 and D e f is some effective "subd- 
iffusion" rate. While the above probability argument gives \i — 1 it is quite possible that 
sticking in the region with k> l\ will give a faster decrease of D + with a higher value of fi. 
As it was discussed in || the growth of pair size should also give logarithmic corrections to 
the coherent localization length in the quasi-ld case ( |1 . 1|) (l c ~ Zi 2 /hr"Zi). 

Another confirmation for the derealization transition below one-particle threshold is 
given by the analysis of the number of effectively excited states. Indeed for v = 3 one 
has AN « <j\<j 1 J_ 2 . According to our data, for sufficiently strong interaction U the quantity 
W = AN/t grows approximately linearly with time (see Fig. 12) while for small U values W 
decreases with time. Contrary to the case v = 2 this indicates that TIP derealization takes 
place for U values bigger than a critical U cr ~ 0.7. Above this critical value the number of 
excited states at a given time grows when k approaches to one-particle derealization border 

k cr (Fig. 13). 



IV. THE EFFECTIVE LLOYD MODEL 

We also studied the model Ql.2|) with the kick perturbation given by Q2.3|) . In this case, 
the non interacting problem can be reduced to the Lloyd model with pseudo-random sites 



energies |fL2j| . For v — 2 and E — the one particle derealization border is k cr w 0.46 |L2 
For TIP problem the behaviour of this model is similar to that of section II. The strong 
enhancement of propagation is demonstrated in Fig. 14. Even if the investigation of this 
model is more difficult for the numerical simulations, our data indicate, as it was in section 
II, that for v = 2 the suppression of a + is always present. In the same way we attribute this 
behaviour to the effective two dimensionality of the model d e ff = 2. 



We also analyzed the tangent model ( |1.2| ), ( |2.3j ) for the case of three frequencies v = 3 



(V(9, 6>i, 02, 3 ) = —2 tan l(2/c(cos6 l + cos^i + cos 62 + cos #3))). This case is similar to 
that discussed in section III. The behaviour of a± is presented in Fig. 15 and indicates the 
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existence of delocalization transition for TIP below one particle derealization border with 
^cv ^ 0.22. 

V. CONCLUSIONS AND DISCUSSIONS 

Our numerical investigations definitely demonstrate the effect of enhancement of the 
localization length for TIP in a random potential. These results were obtained for kicked 
rotators models with frequency modulation. Such approach allows to model efficiently TIP 
problem in an effective dimension d e ff > 2. Numerical data for these models confirm the 
theoretical expectations [f§-|6|] that TIP delocalization in d > 2 is possible below one-particle 
delocalization border. In agreement with || we found TIP pair delocalization and, at the 
same time, a logarithmic growth of the pair size. We attribute this growth to the noise 
produced by the random potential. Indeed a pair propagating in a random potential sees 
different realizations of disorder which act like some effective noise. Such noise originates 
transitions which increase the distance between the two particles. Even if the amplitude 
of these transitions is exponentially decreasing with the two-particle distance, it gives rise 
to logarithmic growth of pair size with time. This in turn produces a subdiffusive pair 
propagation (An + ) 2 ~ D e ft/\ia^t. We give arguments for /i = 1, but it is possible that due 
to sticking in the region with large distance between particles one can have fi > 1. Further 
investigations should be done in order to determine the exact value of /i. Another qualitative 
argument for the subdiffusive propagation is the tunneling between states in which the two 
particles are far from each other (at a distance R ^> l\) and states in which particles stay 
within li (R < li). These states are quasi-degenerate since the spectrum of delocalized 
states is embedded in the spectrum of localized ones. In this situation even an exponentially 
small overlapping between these two kinds of states becomes important and it can lead to 
a subdiffusive pair propagation. Further work should be done for a better understanding of 
the final spectrum structure and the eigenf unctions properties for TIP in 3d. 
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FIGURES 

FIG. 1. One-particle Anderson transition in the model (|2.1| ) with V from ( |2.5| ) v = 2, e = 0.75 
as a function of hopping k. Critical point k cr ~ 1.8. One particle localization lengths (h) and 
diffusion coefficients (-D) are evaluated from the fitting of probability distribution. Error bars 
indicate the standard deviation obtained from an ensemble of 100 (localized) and 10 (diffusive) 
different random realizations. Lines are drawn to fit an eye. 



FIG. 2. Dependence of second moments on time in 2-particles model ( |L2j ) with V from (2.5) 



and v = 2, k = 0.9, e = 0.75,; upper curve is a+(U = 2), middle is <7_ (U = 2), lower is a + (U = 0). 
At t = both particles are at n = n' = 0, basis is —250 < n,n' < 250. Insert shows in larger scale 
the two lower curves in the interval 8 x 10 5 < t < 10 6 . 

FIG. 3. Probability distribution at t = 10 6 as a function of n± = 2 -1 / 2 (n ± n') for the case of 
Fig. 2 with U = 2: P + (n + ) (full line); P_(n_) (dashed); dotted line is the distribution P + (n + ) for 
[7 = 0. 

FIG. 4. Same as Fig.2 but for U = 1. 

FIG. 5. Same as Fig.3 but for U = 1. 

FIG. 6. Dependence of W = AN/t = o+^ct- 1 / 2 ft on time for v = 2, k = 0.9, e = 0.75 and 
U = 2 (full line), U = 1 (dashed line), [7 = 0.5 (dotted line). 

FIG. 7. Dependence of <r+ on A; for v = 2 at f = 2 x 10 5 , [7 = 2 (full circles), U = (open 
circles). 

FIG. 8. Dependence of cr- 1 / 2 on time lnt for cases of Fig.2 (U = 2, k = 0.9) (full upper line) 
and Fig. 4 (U = 1, k = 0.9) (dashed lower line). 
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FIG. 9. One-particle Anderson transition in the model ( |2,1[) with V from ( |3.2|) v = 3, e = 0.9 
as a function of hopping fc. Critical point k cr ~ 1.15. One particle localization lengths (Zi) 
and diffusion coefficient (D) are evaluated from the fitting of probability distribution. Error bars 
indicate the standard deviation for ensemble of 100 (localized) and 10 (diffusive) different random 
realizations. Lines are drawn to fit an eye. 



FIG. 10. Dependence of second moments on time in model ( |1.2| ) with V from ( [3.2| ) and k = 0.7, 
e = 0.9; upper curve is cr+ ( U = 2 ), middle is <r_ (U = 2), lower is <r + ( U = 0). At t = both 
particles are at n = n' = 0, basis is —250 < n,n' < 250, i/ = 3. Inset shows the dependence of 
a- 1 / 2 on ln(t). 

FIG. 11. Probability distribution at t = 1.6 x 10 6 as a function of n± = 2~ l l 2 (n ± n') for the 
case of Fig. 10: P+(n + ) (full line); P_(n_) (dashed); dotted line is the distribution P + (n + ) for 
U = 0. 

FIG. 12. Dependence of W = AN/t = a\ a X J 2 /t on time for v = 3, k = 0.7, e = 0.9 and ?7 = 2 
(upper curve), £/ = 1 (middle curve), U = 0.5 (lower curve). 

FIG. 13. Dependence of <r+ on k for v = 3 at t = 2 x 10 5 , U = 2 (full circles), 17 = (open 
circles). 

FIG. 14. Dependence of second moments on time for the Lloyd model with v = 2, k = 0.35, 
U = 2; upper curve is <r+, lower is <r_. At t = both particles are at n = n' = 0. Basis is 
—256 < n,n' < 256. For {7 = one has a + « 16. 

FIG. 15. Dependence of second moments on time for the Lloyd model with v = 3, k = 0.2, 
{7 = 2; upper curve is cr + , lower is cr_. At i = both particles are at n = n' = 0. Basis is 
128 < n, n! < 128. For {7 = one has a + fs 2. 



15 



2000 
1500 
1000 






n 



6 




6000 
4000 




2000 

<J± 

1000 




